clear;
q=1.6e-19;
m=1.6725e-27;
V0=1e6;sita=pi/4;
Vy=0,Vz=V0*cos(sita),Vx=V0*sin(sita),
c=[0,1.1,0,Vx,Vy,Vz];%给定初始量
[t,y]=ode45('cjxxfun',[0:1e-7:1.5e-6],c,[],q,m);%解微分方程组
plot3(y(:,1),y(:,2),y(:,3),'linewidth',2);%绘制三维粒子运动轨迹
xlabel('z');ylabel('y');zlabel('x');
grid on;title('磁镜现象')
figure
subplot(2,1,1);
plot(y(:,1),y(:,2),'linewidth',2);%从yoz面看到的轨迹
xlabel('z');ylabel('y');
grid on;title('yoz平面')
subplot(2,1,2);
plot(y(:,1),y(:,3),'linewidth',2);%从xoz面观察的轨迹
xlabel('z');ylabel('x');
title('xoz平面');grid on;


